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Breit— Wigner Enhancement of Dark Matter Annihilation 
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We point out that annihilation of dark matter in the galactic halo can be enhanced relative to 
that in the early universe due to a Breit-Wigner tail, if the dark matter annihilates through a pole 
just below the threshold. This provides a new explanation to the "boost factor" which is suggested 
by the recent data of the PAMELA, ATIC and PPB-BETS cosmic-ray experiments. 



Dark matter of the universe has been discussed since 
1933, yet its nature still remains elusive Seventy-five 
years later we only managed to restrict its mass between 
10"'^^ and 10^" GeV, demonstrating our lack of under- 
standing. However, the thermal relic of a Weakly Inter- 
acting Massive Particle (WIMP) is arguably best theo- 
retically motivated because it has the same mass scale 
as the anticipated new physics that would explain why 
our universe is a superconductor (electroweak symme- 
try breaking) . Hopes are high to discover WIMPs at the 
forthcoming LHC experiments, to detect them directly in 
sensitive underground experiments, as well as to observe 
signals of WIMP annihilations from the galactic center 
or the halo in high-energy cosmic rays. 

Recent observations of the PAMELA 01, ATIC and 
PPB-BETS [4] experiments strongly suggest the exis- 
tence of a new source of positrons (and electrons) in cos- 
mic rays. The most interesting interpretation of these 
results is the annihilation of the dark matter with a mass 
at the TeV scale. However, such interpretation requires 
that the annihilation cross section of the dark matter in 
the galactic halo is much larger (by a factor of O(IOO)) 
than the one appropriate to explain the dark matter relic 
density precisely measured by the WMAP experiment 'B'l . 

The enhancement of the dark matter annihilation in 
the galactic halo is called a "boost factor." So far, there 
have been several proposals to explain the origin of the 
boost factor both from astrophysics such as the enhanced 
local dark matter density, and from particle physics such 
as the Sommerfeld enhancement due to an attractive 
force among the dark matter particles ^6;] . 

In this letter, we propose a new explanation of the 
boost factor. We consider the dark matter which annihi- 
lates via a narrow Breit-Wigner resonance just below the 
threshold. When the resonance mass is just below twice 
the dark matter mass, the annihilation cross section be- 
comes sensitive to the velocity of the dark matter. In such 
a case, the time evolution of the dark matter abundance 
is quite different from the one in the usual non-resonant 
case of annihilation, and we find that the annihilation 
cross section in the halo is enhanced compared to what 
is inferred from the relic abundance. As we will show, 



the cross section required from the dark matter density 
can be large enough to explain the PAMELA, ATIC, and 
PPB-BETS results, and hence, we do not need in our 
proposal any additional boost factor due to an overdense 
region in the halo or the Sommerfeld enhancement. 

Cross Section Just Above a Pole. In this study, we 
assume that the dark matter with mass m annihilates 
via a narrow resonance R. For a simplicity, we consider 
a scalar resonance, although generalization to arbitrary 
spins is straight forward. The general formula for the 
scattering cross section via a resonance R is given by 

167r M^F^ 

where M and F are the mass and the decay rate of the 
resonance i?, respectively. Two body initial and final 
states are assumed and Pi = -^/l — Am? /M"^ is the initial 
state phase space factor evaluated on the resonance while 
fii = a/1 — im'^/E^^ at the center of mass energy of 
the collision. The branching fractions of the resonance 
into the initial and final states are denoted by Bi and 
Bf, respectively. Note that since we are assuming an 
unphysical pole, i.e., 2m > M, Bi and f3i are not physical 
and should be understood as analytic continuations of 
those quantities from the physical region, 2m < M . Even 
so, a combination Bi/ Pi is well-defined in both regions, 
and hence, the above cross section is also well-defined 
even in the unphysical pole case. 

The dark matter annihilation in the early universe 
must be thermally averaged. On the other hand, the 
dark matter annihilation in the galactic halo is averaged 
over the velocity distribution, which can be approximated 
by the Maxwellian distribution. In both cases, the dark 
matter is non-relativistic. Therefore, in either case, we 
can use the Gaussian average, 

{(Tv,,y) - — ^ /dz/id«2e-(""+"")/''^"2aA, (2) 



where vi^2 are the velocities of the initial states and we 
have used Wrci = 2/?.^. The non-relativistic approximation 
gives 



EI^^ = 4m^ + m^{v^c\f, {v^c\ = vi - V2). 



(3) 
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Now, let us consider the annihilation process near a 
narrow resonance, i.e., 

m2 = W(1-(5), |(5| < 1. (4) 

Note that positive S implies that the pole is just below 
the threshold of the dark matter annihilation. With this 
notation, we can rewrite the above cross section as, 
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where we have defined, 

7 = T/M. 
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(5) 
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Furthermore, we have verified that we can approximate 
the Gaussian integral reasonably well by 
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M2^, ((5 + ew2)2 +72 
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(7) 



where a parameter ^ « l/"\/2 gives the best fit to the 
numerical results for <SC 1 and S > 0. This expression 
shows that the denominator is dominated by the v^^i term 
for I (5 1, 7 <C ^Vq, while the other terms dominate when 
the velocity is much smaller. Therefore, the cross sec- 
tion is sensitive to the parameters S and 7, and have an 
enhanced behavior at the lower temperature for small 6 
and 7. This main point of this letter can be seen easily 
in a schematic plot in Fig. [TJ 

Note that the approximation given in Eq. ([7]) is not a 
good one for the dark matter with a rather large veloc- 
ity, i.e., vo ^ 0.1 Q- The approximation also becomes 
worse around the pole in the physical region, possible if 
5 < 0. In the following analysis, we mainly consider the 
unphysical pole (5 > 0. As we will see below, for this case, 
the resultant dark matter density is mostly determined 
by the dynamics of the dark matter at <C 1, and hence 
the approximated cross section works quite well. We will 
later also briefly discuss the case of a physical pole, S < 0. 

Time Evolution of Dark Matter Density. Now let us 
consider the time evolution of the relic density of the dark 
matter. The most notable feature of our proposal is that 
the annihilation process does not freeze out even after 
the "freeze-out" time for the usual non-resonant annihi- 
lation models. That is, the interaction rate can be larger 
than the Hubble expansion rate even after a "freeze-out" 
time, since the cross section is enhanced at the lower 
temperature. In the following, to avoid confusion, we 
define the "freeze-out" time 55/ by the usual meaning at 
which the yield of the dark matter Y — n/s starts de- 
viating from the value in the thermal equilibrium, i.e., 
Y — Yeq — 0{Yeq). As we will show, the actual freeze- 
out time is much later than i/ . In this section, we mainly 
consider the unphysical pole, S > 0. 

Following Ref. [11, section 5.2, we write down the 
Boltzmann equation of the yield of the dark matter [lo| , 

' "'■'■•'^ (Y'-Yi,, (8) 



dY 
dx 
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FIG. 1: A schematic plot that shows dispersion in relative 
velocity v^^i and an unphysical pole in the cross section at 
Wrei < (below threshold) . It is clear that a smaller dispersion 
Wo gives a larger overlap with the Breit-Wigner tail in the 
cross section and hence an enhanced averaged cross section. 



Here, we have used the following definitions 
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The parameter (gi) is the number of the degrees 
of freedom for massless particles (dark matter), respec- 
tively. Note that we have used the reduced Planck scale 
MpL ~ 2.4 X 10^8 GeV. 

As we have defined, the "freeze-out" time i/ is deter- 
mined by y — Ieq — 0(1eq), and hence, the value of i/ 
is not so sensitive to the parameters S and 7 and mainly 
determined by the exponential suppression factor in Ieq- 
Thus, we can expect that the value oi x f is comparable to 
the freeze-out temperature Xf in the usual non-resonant 
annihilation models, i.e., xj ^ Xf ~ 0(10). 

Unlike the non-rcsonant case where the annihilation 
cross section stays constant (crwroi) once non-relativistic, 
however, the annihilation cross section here increases as 
temperature drops. As a result, the annihilation process 
does not freeze out even for x > xj and the dark matter 
keeps annihilating until the temperature comes down to 



Xb 



^ "'^ X max [(5, 7] ^ ^ Xf. 
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Below this temperature, Eq. 

dY _ A 
dx 



reduces to 



r2 



(10) 



(11) 
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and we obtain an asymptotic solution, 

Yoo — \ Xb ^ \ X max [5, -y]"^ . 

A A 



(12) 



In the Fig.[2l we show the time evolution of the yield 
Y of the dark matter for a given parameter set. Here, 
we have used the numerical result of the Gaussian av- 
erage of the cross section We also show the time 



evolution of the yield with the approximate cross sec- 
tion given in Eq. d?]), Vapp- As we have expected, the 
yield deviates from Yeq at if = 0(10), while the actual 
freeze-out occurs at Xf, ^ if. From the figure, we see 
that Y < Yapp during if < x < x^- This means that 
the averaged cross section at 5/ is somewhat larger than 
the approximate one, while the final result is determined 
by the late time dynamics where the approximate cross 
section works well. The figure shows that the above ap- 
proximate asymptotic solution gives a good estimate of 
the resultant yield of the dark matter. 

It is worth comparing the asymptotic solution in 
Eq. (|12p with the asymptotic solution in the usual non- 
resonant (5- wave) annihilation models [sl, 



Y^ 



1 



(13) 



with Xf = 0(10). These two solutions show that the dark 
matter abundance in our proposal is larger by a factor 



Xb_ 



max [5, 7] 
0(10) 



(14) 



when we assume the same cross section at the zero tem- 
perature for both models. In the Fig.[5J we have showed 
the time evolution of the yield in the usual non-resonant 
annihilation models assuming the same cross section at 
the zero temperature (Tq used in Y {i.e., A — 10^). As 
expected, the yield in the non-resonant annihilation is 
more suppressed compared to the resonant annihilation. 
The physical reasons of this enhancement are as follows. 
First, in our proposal, the cross section is suppressed by 
{if /xbY at if, which results in a relatively larger abun- 
dance in a period of x > if. Second, the annihilation 
process is relatively less effective during if < x < x^ 
compared to the usual annihilation case, although the 
annihilation does not freeze out in that period. 

Before closing this section, we mention the model with 
a physical pole, i.e., 5 < 0. In this case, the cross sec- 
tion given in Eq. ([7]) poorly approximates the thermal av- 
eraged cross section [7| . Especially, the thermal average 
can pick up the cross section at the pole, v'^^^ = A\d\ when 
the temperature is rather high, i.e., x~~^ ^ \S\. Thus, the 
averaged cross section can be much higher than that ex- 
pected in the unphysical pole where the cross section is 
suppressed by x'^{6'^+j'^) at high temperature. Thus, the 
annihilation cross section at if , which typically satisfies 
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FIG. 2: The time evolution of the yield Y of the dark mat- 
ter in terms of the parameter x = m/T for given values of 5 
and 7 (the solid line). In this figure, we assume 5 > and 
the pole is not in the physical region. The long-dashed line 
labeled Kapp is the evolution with the approximated cross sec- 
tion in Eq. ([7]). The dashed line labeled Yoo is the asymptotic 
solution Yoc given in Eq. ((12} . The short-dashed line repre- 
sents the equilibrium yield Yeq- The dash-dotted line labeled 
ynon-ros shows the time evolution of the yield in the usual 
(non-resonant) freeze-out assuming the same cross section at 
the low temperature ctq used in Y (see Eq. (|13p ). The boost 
factor is the asymptotic value of K/ynon-rcs. 
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FIG. 3; The time evolution of the yield of the dark matter Y 
for S < (the solid line). Everything else is the same as in 
Fig. [2] There is practically no boost for this parameter set. 



ij^ 3> |<^|,7, is much larger than that in the case of the 
unphysical pole. Therefore, the yield enhancement seen 
in the unphysical pole is much smaller in this case. 

Effective Boost Factor. Finally, let us work out the 
boost factor for annihilation in the galactic halo. As 
we have seen, the dark matter abundance is enhanced 
by a factor of {xb/xf) for 5 > compared to the usual 
non-resonant annihilation models for a given annihilation 
cross section at the zero temperature, uo- This means 
that the cross section ctq must be larger than the one 
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expected in the usual models to reproduce the observed 
dark matter density. 

More explicitly, the yield given in Eq. is translated 
to the mass density parameter, 



0.1 X 



where we have used xj 
the observed dark matter density 17 /i^ 



10-9 GeV 

CTO 

~ 20 and 



(15) 



CTo ~ 10^9 X 



200. Therefore, 
i 0.1 requires 

(16) 



which is much larger than what is expected in the usual 
annihilation case, i.e., 



CTo ~ 10-9 GeV" 



(17) 



In the galactic halo, the average velocity is given 
by vq ~ IQ-'^, and the cross section is well approxi- 
mated by the one at the zero temperature as long as 
Vq <C (5, 7. Thus, we can achieve the large annihilation 
cross section suggested by the PAMELA, PPB-BETS, 
and ATIC anomalies, (CTWrci)T=o = 0(10-^^-''^) GcY-^ 
for (5, 7 < 10"^. Therefore, in our proposal, we can ex- 
plain the large annihilation cross section in the galactic 
halo without other boost factors. 

For convenience, we could define an effective boost fac- 
tor as the ratio between the cross sections in the usual 
and our models. From Eqs. (fT6|) and ((T7)) we obtain 



BF^^ 



max [6, 7] 



0(10) 



(18) 



In this way, we can explain the boost factor of O(IOO), 
in a model with S,j — O{10~^), as seen in Fig. 21 

Discussion. In this letter, we showed that the boost 
factor required in recent observations of cosmic ray elec- 
trons and positrons can be obtained if the dark matter 
annihilates via a narrow resonance just below the thresh- 
old. Because the annihilation in the early universe is 
suppressed by the Breit-Wigner tail, the observed dark 
matter density requires a larger-than-normal cross sec- 
tion which can be consistent with the PAMELA, ATIC 
and PPB-BETS results. The required cross section is 
achieved for im"^ - M^, MT < lO^^lVP. 

Note that the small decay width F <C lO^^M can be 
achieved rather easily, when the resonance i? is a weakly 
coupled particle. In some models, a coincidence in masses 
can also be naturally realized. For example, if the dark 
matter is the lightest Kaluza-Klein (KK) particle and the 
resonance is at the second KK level, we have the relation 
M = 2m at the tree-level [i] [i3 (see also a recent dis- 
cussion in Ref. [1]). Thus, in such models, we may well 
have a small mass splitting 5 as a result of small radiative 
corrections to the tree-level mass relation. 
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FIG. 4; The boost factor on the (5, 7) plane. Thermal aver- 
age is done numerically without relying on the approximation 
Eq. (0. 
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